Design and synthesis of Fsp3-enriched spirocyclic-substituted diarylpyrimidine derivatives as novel HIV-1 NNRTIs.

In this study, a novel series of diarylpyrimidine derivatives with Fsp3-enriched spirocycles were designed and synthesized to further explore the chemical space of the hydrophobic channel of the NNRTI-binding pocket. The biological evaluation results showed that most of the compounds displayed effective inhibitory potency against the HIV-1 wild-type strain, with EC50 values ranging from micromolar to submicromolar levels. Among them, TT6 turned out to be the most effective inhibitor with an EC50 value of 0.17 µM, demonstrating up to 47 times more active than that of reference drug 3TC (EC50 = 8.01 µM). More encouragingly, TT6 was found to potently inhibit the HIV-1 mutant strain K103N with an EC50 value of 0.69 µM, being about 6-fold more potent than 3TC (EC50 = 3.68 µM) and NVP (EC50 = 4.62 µM). Furthermore, TT6 exhibited the most potent inhibitory activity toward HIV-1 reverse transcriptase with an IC50 value of 0.33 µM. Additionally, molecular simulation studies were conducted to investigate the binding modes between TT6 and NNRTI-binding pocket, which may provide valuable clues for the follow-up structural optimizations.

Selected Milestones in Antiviral Drug Development.

This review article will describe the (wide) variety of approaches that I envisaged to develop a specific therapy for viral infections: (i) interferon and its inducers, (ii) HSV, VZV and CMV inhibitors, (iii) NRTIs (nucleoside reverse transcriptase inhibitors), NtRTIs (nucleotide reverse transcriptase inhibitors) and NNRTIs (non-nucleoside reverse transcriptase inhibitors) as HIV inhibitors, (iv) NtRTIs as HBV inhibitors, and finally, (v) the transition of an HIV inhibitor to a stem cell mobilizer, as exemplified by AMD-3100 (Mozobil®).

HIV-1 Reverse Transcriptase Expression in HPV16-Infected Epidermoid Carcinoma Cells Alters E6 Expression and Cellular Metabolism, and Induces a Hybrid Epithelial/Mesenchymal Cell Phenotype.

The high incidence of epithelial malignancies in HIV-1 infected individuals is associated with co-infection with oncogenic viruses, such as high-risk human papillomaviruses (HR HPVs), mostly HPV16. The molecular mechanisms underlying the HIV-1-associated increase in epithelial malignancies are not fully understood. A collaboration between HIV-1 and HR HPVs in the malignant transformation of epithelial cells has long been anticipated. Here, we delineated the effects of HIV-1 reverse transcriptase on the in vitro and in vivo properties of HPV16-infected cervical cancer cells. A human cervical carcinoma cell line infected with HPV16 (Ca Ski) was made to express HIV-1 reverse transcriptase (RT) by lentiviral transduction. The levels of the mRNA of the E6 isoforms and of the factors characteristic to the epithelial/mesenchymal transition were assessed by real-time RT-PCR. The parameters of glycolysis and mitochondrial respiration were determined using Seahorse technology. RT expressing Ca Ski subclones were assessed for the capacity to form tumors in nude mice. RT expression increased the expression of the E6*I isoform, modulated the expression of E-CADHERIN and VIMENTIN, indicating the presence of a hybrid epithelial/mesenchymal phenotype, enhanced glycolysis, and inhibited mitochondrial respiration. In addition, the expression of RT induced phenotypic alterations impacting cell motility, clonogenic activity, and the capacity of Ca Ski cells to form tumors in nude mice. These findings suggest that HIV-RT, a multifunctional protein, affects HPV16-induced oncogenesis, which is achieved through modulation of the expression of the E6 oncoprotein. These results highlight a complex interplay between HIV antigens and HPV oncoproteins potentiating the malignant transformation of epithelial cells.

Synthesis and anti-HIV activities of phorbol derivatives.

In this study, 37 derivatives of phorbol esters were synthesized and their anti-HIV-1 activities evaluated, building upon our previous synthesis of 51 phorbol derivatives. 12-Para-electron-acceptor-trans-cinnamoyl-13-decanoyl phorbol derivatives stood out, demonstrating remarkable anti-HIV-1 activities and inhibitory effects on syncytia formation. These derivatives exhibited a higher safety index compared with the positive control drug. Among them, 12-(trans-4-fluorocinnamoyl)-13-decanoyl phorbol, designated as compound 3c, exhibited the most potent anti-HIV-1 activity (EC50 2.9 nmol·L-1, CC50/EC50 11 117.24) and significantly inhibited the formation of syncytium (EC50 7.0 nmol·L-1, CC50/EC50 4891.43). Moreover, compound 3c is hypothesized to act both as an HIV-1 entry inhibitor and as an HIV-1 reverse transcriptase inhibitor. Isothermal titration calorimetry and molecular docking studies indicated that compound 3c may also function as a natural activator of protein kinase C (PKC). Therefore, compound 3c emerges as a potential candidate for developing new anti-HIV drugs.

Proposal of pharmacophore model for HIV reverse transcriptase inhibitors: Combined mutational effect analysis, molecular dynamics, molecular docking and pharmacophore modeling study.

OBJECTIVES: Antiretroviral therapy (ART) efficacy is jeopardized by the emergence of drug resistance mutations in HIV, compromising treatment effectiveness. This study aims to propose novel analogs of Effavirenz (EFV) as potential direct inhibitors of HIV reverse transcriptase, employing computer-aided drug design methodologies. METHODS: Three key approaches were applied: a mutational profile study, molecular dynamics simulations, and pharmacophore development. The impact of mutations on the stability, flexibility, function, and affinity of target proteins, especially those associated with NRTI, was assessed. Molecular dynamics analysis identified G190E as a mutation significantly altering protein properties, potentially leading to therapeutic failure. Comparative analysis revealed that among six first-line antiretroviral drugs, EFV exhibited notably low affinity with viral reverse transcriptase, further reduced by the G190E mutation. Subsequently, a search for EFV-similar inhibitors yielded 12 promising molecules based on their affinity, forming the basis for generating a pharmacophore model. RESULTS: Mutational analysis pinpointed G190E as a crucial mutation impacting protein properties, potentially undermining therapeutic efficacy. EFV demonstrated diminished affinity with viral reverse transcriptase, exacerbated by the G190E mutation. The search for EFV analogs identified 12 high-affinity molecules, culminating in a pharmacophore model elucidating key structural features crucial for potent inhibition. CONCLUSION: This study underscores the significance of EFV analogs as potential inhibitors of HIV reverse transcriptase. The findings highlight the impact of mutations on drug efficacy, particularly the detrimental effect of G190E. The generated pharmacophore model serves as a pivotal reference for future drug development efforts targeting HIV, providing essential structural insights for the design of potent inhibitors based on EFV analogs identified in vitro.

Internal RNA 2'-O-methylation on the HIV-1 genome impairs reverse transcription.

Viral RNA genomes are modified by epitranscriptomic marks, including 2'-O-methylation that is added by cellular or viral methyltransferases. 2'-O-Methylation modulates RNA structure, function and discrimination between self- and non-self-RNA by innate immune sensors such as RIG-I-like receptors. This is illustrated by human immunodeficiency virus type-1 (HIV-1) that decorates its RNA genome through hijacking the cellular FTSJ3 2'-O-methyltransferase, thereby limiting immune sensing and interferon production. However, the impact of such an RNA modification during viral genome replication is poorly understood. Here we show by performing endogenous reverse transcription on methylated or hypomethylated HIV-1 particles, that 2'-O-methylation negatively affects HIV-1 reverse transcriptase activity. Biochemical assays confirm that RNA 2'-O-methylation impedes reverse transcriptase activity, especially at low dNTP concentrations reflecting those in quiescent cells, by reducing nucleotide incorporation efficiency and impairing translocation. Mutagenesis highlights K70 as a critical amino acid for the reverse transcriptase to bypass 2'-O-methylation. Hence, the observed antiviral effect due to viral RNA 2'-O-methylation antagonizes the FTSJ3-mediated proviral effects, suggesting the fine-tuning of RNA methylation during viral replication.

Current insights and molecular docking studies of HIV-1 reverse transcriptase inhibitors.

Human immunodeficiency virus (HIV) causes acquired immunodeficiency syndrome (AIDS), a lethal disease that is prevalent worldwide. According to the Joint United Nations Programme on HIV/AIDS (UNAIDS) data, 38.4 million people worldwide were living with HIV in 2021. Viral reverse transcriptase (RT) is an excellent target for drug intervention. Nucleoside reverse transcriptase inhibitors (NRTIs) were the first class of approved antiretroviral drugs. Later, a new type of non-nucleoside reverse transcriptase inhibitors (NNRTIs) were approved as anti-HIV drugs. Zidovudine, didanosine, and stavudine are FDA-approved NRTIs, while nevirapine, efavirenz, and delavirdine are FDA-approved NNRTIs. Several agents are in clinical trials, including apricitabine, racivir, elvucitabine, doravirine, dapivirine, and elsulfavirine. This review addresses HIV-1 structure, replication cycle, reverse transcription, and HIV drug targets. This study focuses on NRTIs and NNRTIs, their binding sites, mechanisms of action, FDA-approved drugs and drugs in clinical trials, their resistance and adverse effects, their molecular docking studies, and highly active antiretroviral therapy (HAART).

Variable antiviral activity of islatravir against M184I/V mutant HIV-1 selected during antiretroviral therapy.

BACKGROUND: Islatravir is a new antiretroviral drug that inhibits the reverse transcriptase (RT) of HIV-1 through multiple mechanisms. It is proposed to be used in combination with doravirine, a new NNRTI. M184V/I mutations have been shown to reduce the in vitro antiviral activity of islatravir, but their effect when pre-selected during ART has not been investigated. METHODS: HIV-1 rt sequences were obtained from four individuals of the Garrahan HIV cohort prior to, or during virological failure to ART. HIV-1 infectious molecular clones were constructed on an NL4-3 backbone, and infectious viruses were produced by transfection of 293T cells. Fold-changes in IC50 were calculated for each mutant versus the NL4-3 WT. HIV-1 phenotypic drug resistance was tested in vitro against NRTIs and NNRTIs. RESULTS: In all the cases, M184I/V, either alone or in the presence of other mutations, was associated with reduced susceptibility to islatravir, abacavir and lamivudine. Viruses carrying M184V/I showed variable levels of resistance to islatravir (4.8 to 33.8-fold). The greatest reduction in susceptibility was observed for viruses carrying the mutations M184V + V106I (33.8-fold resistance) or M184V + I142V (25.2-fold resistance). For NNRTIs, the presence of V106I alone did not affect susceptibility to doravirine or etravirine, but showed a modest reduction in susceptibility to efavirenz (6-fold). Susceptibility to doravirine was slightly reduced only for one of the mutants carrying V106I in combination with Y181C and M184V. CONCLUSIONS: Mutations and polymorphisms selected in vivo together with M184V/I depend on the viral genetic context and on ART history, and could affect the efficacy of islatravir once available for use in the clinic.

Targeting the Structural Maturation Pathway of HIV-1 Reverse Transcriptase.

Formation of active HIV-1 reverse transcriptase (RT) proceeds via a structural maturation process that involves subdomain rearrangements and formation of an asymmetric p66/p66' homodimer. These studies were undertaken to evaluate whether the information about this maturation process can be used to identify small molecule ligands that retard or interfere with the steps involved. We utilized the isolated polymerase domain, p51, rather than p66, since the initial subdomain rearrangements are largely limited to this domain. Target sites at subdomain interfaces were identified and computational analysis used to obtain an initial set of ligands for screening. Chromatographic evaluations of the p51 homodimer/monomer ratio support the feasibility of this approach. Ligands that bind near the interfaces and a ligand that binds directly to a region of the fingers subdomain involved in subunit interface formation were identified, and the interactions were further characterized by NMR spectroscopy and X-ray crystallography. Although these ligands were found to reduce dimer formation, further efforts will be required to obtain ligands with higher binding affinity. In contrast with previous ligand identification studies performed on the RT heterodimer, subunit interface surfaces are solvent-accessible in the p51 and p66 monomers, making these constructs preferable for identification of ligands that directly interfere with dimerization.

The Application of Prodrugs as a Tool to Enhance the Properties of Nucleoside Reverse Transcriptase Inhibitors.

Antiretroviral Therapy (ART) is an effective treatment for human immunodeficiency virus (HIV) which has transformed the highly lethal disease, acquired immunodeficiency syndrome (AIDS), into a chronic and manageable condition. However, better methods need to be developed for enhancing patient access and adherence to therapy and for improving treatment in the long term to reduce adverse effects. From the perspective of drug discovery, one promising strategy is the development of anti-HIV prodrugs. This approach aims to enhance the efficacy and safety of treatment, promoting the development of more appropriate and convenient systems for patients. In this review, we discussed the use of the prodrug approach for HIV antiviral agents and emphasized nucleoside reverse transcriptase inhibitors. We comprehensively described various strategies that are used to enhance factors such as water solubility, bioavailability, pharmacokinetic parameters, permeability across biological membranes, chemical stability, drug delivery to specific sites/organs, and tolerability. These strategies might help researchers conduct better studies in this field. We also reported successful examples from the primary therapeutic classes while discussing the advantages and limitations. In this review, we highlighted the key trends in the application of the prodrug approach for treating HIV/AIDS.

Web Service for HIV Drug Resistance Prediction Based on Analysis of Amino Acid Substitutions in Main Drug Targets.

Predicting viral drug resistance is a significant medical concern. The importance of this problem stimulates the continuous development of experimental and new computational approaches. The use of computational approaches allows researchers to increase therapy effectiveness and reduce the time and expenses involved when the prescribed antiretroviral therapy is ineffective in the treatment of infection caused by the human immunodeficiency virus type 1 (HIV-1). We propose two machine learning methods and the appropriate models for predicting HIV drug resistance related to amino acid substitutions in HIV targets: (i) k-mers utilizing the random forest and the support vector machine algorithms of the scikit-learn library, and (ii) multi-n-grams using the Bayesian approach implemented in MultiPASSR software. Both multi-n-grams and k-mers were computed based on the amino acid sequences of HIV enzymes: reverse transcriptase and protease. The performance of the models was estimated by five-fold cross-validation. The resulting classification models have a relatively high reliability (minimum accuracy for the drugs is 0.82, maximum: 0.94) and were used to create a web application, HVR (HIV drug Resistance), for the prediction of HIV drug resistance to protease inhibitors and nucleoside and non-nucleoside reverse transcriptase inhibitors based on the analysis of the amino acid sequences of the appropriate HIV proteins from clinical samples.

Biochemical and structural comparisons of non-nucleoside reverse transcriptase inhibitors against feline and human immunodeficiency viruses.

BACKGROUND: Feline immunodeficiency virus (FIV) causes an acquired immunodeficiency-like syndrome in cats. FIV is latent. No effective treatment has been developed for treatment the infected cats. The first and second generations non-nucleoside reverse transcriptase inhibitors (NNRTIs) for HIV treatment, nevirapine (NVP) and efavirenz (EFV), and rilpivirine (RPV), were used to investigate the potential of NNRTIs for treatment of FIV infection. OBJECTIVE: This study aims to use experimental and in silico approaches to investigate the potential of NNRTIs, NVP, EFV, and RPV, for inhibition of FIV reverse transcriptase (FIV-RT). METHODS: The FIV-RT and human immunodeficiency virus reverse transcriptase (HIV-RT) were expressed and purified using chromatography approaches. The purified proteins were used to determine the IC50 values with NVP, EFV, and RPV. Surface plasmon resonance (SPR) analysis was used to calculate the binding affinities of NNRTIs to HIV-RT and FIV-RT. The molecular docking and molecular dynamic simulations were used to demonstrate the mechanism of FIV-RT and HIV-RT with first and second generation NNRTI complexes. RESULTS: The IC50 values of NNRTIs NVP, EFV, and RPV against FIV-RT were in comparable ranges to HIV-RT. The SPR analysis showed that NVP, EFV, and RPV could bind to both enzymes. Computational calculation also supports that these NNRTIs can bind with both FIV-RT and HIV-RT. CONCLUSIONS: Our results suggest the first and second generation NNRTIs (NVP, EFV, and RPV) could inhibit both FIV-RT and HIV-RT.

In situ click chemistry-based discovery of 1,2,3-triazole-derived diarylpyrimidines as novel HIV-1 NNRTIs by exploiting the tolerant region I in binding pocket.

HIV-1 reverse transcriptase (RT) is considered as one of the most significant targets for the anti-HIV-1 drug design due to their determined mechanism and well-decoded crystal structure. As a part of our continuous efforts towards the development of potent HIV-1 non-nucleoside reverse transcriptase inhibitors (NNRTIs) by exploiting the tolerant region I of NNRTIs binding pocket (NNIBP), the miniaturized parallel synthesis via CuAAC click chemistry reaction followed by in situ biological screening have been performed in this work. The in situ enzyme inhibition screening results showed that 14 compounds exhibited higher or equivalent inhibitory activity compared to the lead K-5a2 and ETR. Anti-HIV-1 activity results indicated that C1N51 displayed the most potent activity (EC50 = 0.01-0.26 µM) against wild-type and a panel of NNRTIs-resistant strains. Moreover, the molecular simulation demonstrated that the newly introduced triazole ring could develop new hydrogen bonds with Lys103 and Pro236, which explained the feasibility of introducing triazole in the tolerant region I of the RT binding pocket.

Exploring novel HIV-1 reverse transcriptase inhibitors with drug-resistant mutants: A double mutant surprise.

HIV-1 reverse transcriptase (RT) remains a key target for HIV drug development. As successful management of the disease requires lifelong treatment, the emergence of resistance mutations is inevitable, making development of new RT inhibitors, which remain effective against resistant variants crucial. To this end, previous computationally guided drug design efforts have resulted in catechol diether compounds, which inhibit wildtype RT with picomolar affinities and appear to be promising preclinical candidates. To confirm that these compounds remain potent against Y181C, a widespread mutation conferring resistance to first generation inhibitors, they were screened against the HIV-1 N119 clinical isolate, reported as a Y181C single mutant. In comparison to a molecular clone with the same mutation, N119 appears less susceptible to inhibition by our preclinical candidate compounds. A more detailed sequencing effort determined that N119 was misidentified and carries V106A in combination with Y181C. While both indolizine and naphthalene substituted catechol diethers are potent against the classical Y181C single mutant, the addition of V106A confers more resistance against the indolizine derivatives than the naphthalene derivatives. Crystal structures presented in this study highlight key features of the naphthyl group, which allow these compounds to remain potent in the double mutant, including stronger interactions with F227 and less reliance on V106 for stabilization of the ethoxy-uracil ring, which makes critical hydrogen bonds with other residues in the binding pocket.

Strategies in the Design and Development of Non-Nucleoside Reverse Transcriptase Inhibitors (NNRTIs).

AIDS (acquired immunodeficiency syndrome) is a potentially life-threatening infectious disease caused by human immunodeficiency virus (HIV). To date, thousands of people have lost their lives annually due to HIV infection, and it continues to be a big public health issue globally. Since the discovery of the first drug, Zidovudine (AZT), a nucleoside reverse transcriptase inhibitor (NRTI), to date, 30 drugs have been approved by the FDA, primarily targeting reverse transcriptase, integrase, and/or protease enzymes. The majority of these drugs target the catalytic and allosteric sites of the HIV enzyme reverse transcriptase. Compared to the NRTI family of drugs, the diverse chemical class of non-nucleoside reverse transcriptase inhibitors (NNRTIs) has special anti-HIV activity with high specificity and low toxicity. However, current clinical usage of NRTI and NNRTI drugs has limited therapeutic value due to their adverse drug reactions and the emergence of multidrug-resistant (MDR) strains. To overcome drug resistance and efficacy issues, combination therapy is widely prescribed for HIV patients. Combination antiretroviral therapy (cART) includes more than one antiretroviral agent targeting two or more enzymes in the life cycle of the virus. Medicinal chemistry researchers apply different optimization strategies including structure- and fragment-based drug design, prodrug approach, scaffold hopping, molecular/fragment hybridization, bioisosterism, high-throughput screening, covalent-binding, targeting highly hydrophobic channel, targeting dual site, and multi-target-directed ligand to identify and develop novel NNRTIs with high antiviral activity against wild-type (WT) and mutant strains. The formulation experts design various delivery systems with single or combination therapies and long-acting regimens of NNRTIs to improve pharmacokinetic profiles and provide sustained therapeutic effects.

The impact of the M184V resistance mutation on treatment outcomes in patients with HIV infection: a systematic review and meta-analysis.

HIV is a global deliberating infectious disease. Of note, more than 36 million people living with HIV (PLHIV) with approximately newly diagnosed 1.5 million cases annually. M184V is a single base mutation in the highly conserved YMDD domain of reverse transcriptase (RT). It is one of the most encountered resistances associated with mutations to nucleoside RT inhibitors. There were continuous efforts to evaluate the impact of M184V mutation on the treatment outcomes in PLHIV. Therefore, the present systematic review was executed to reveal the virological failure, virological suppression, and resistance to antiretroviral therapy (ART) regimens in PLHIV with the M184V mutation. All clinical studies comparing the treatment outcomes among PLHIV harboring or not harboring M184V mutation were appropriate for systematic review and meta-analysis. The present systematic review included six articles, encompassing 4760 PLHIV. Of them, 1222 (25.67%) patients had M184V mutation, while 3538 (74.32%) PLHIV did not. The meta-analysis showed that patients with M184V mutation were 1.87 times more liable to virological failure (risk ratio [RR] 1.87; 95% 1.09, 3.20; p = 0.02). Furthermore, pooling the data from two studies revealed a significantly higher risk of viral blips (RR 2.26; 95% 1.47, 3.46; p = 0.0002). Concerning discontinuation of ART, there was no statistical difference between patients with and without M184V mutation (RR: 0.99; 95% 0.78, 1.25; p = 0.90). The present study revealed the negative impact of the M184V mutation on treatment outcomes in PLHIV. This included a higher risk of virological failure and viral blips, relative to patients without the mutation. Such patients may benefit from more aggressive and combined therapy for better disease management.

4-Phenylcoumarin derivatives as new HIV-1 NNRTIs: Design, synthesis, biological activities, and computational studies.

A series of 4-phenylcoumarin derivatives was synthesized and evaluated for their cellular anti-HIV-1 and HIV-2 activities as well as their inhibitory effects against HIV-1 reverse transcriptase (RT). The hydrazone compound 8b and the ethylthiosemicarbazide derivative 4c showed the best inhibition activity against wild-type (WT) HIV-1. The promising compounds were further evaluated against HIV-1 RT and exhibited significant inhibitory activity with compound 8b showing comparable effect to the reference NNRTI Efavirenz (IC50 = 9.01 nM). Structure activity relationship study revealed the importance of 6-chloro and 4-phenyl substituents for optimum activity, as well as the 5-atoms linker (=N-NH-CO-CH2-O-) at position 7 of coumarin scaffold that can support the rotation and flexibility of compound 8b to fit well in the binding pocket. The molecular docking of compound 8b demonstrated a typical seahorse binding mode with better binding interactions that covered more residues when compared to Efavirenz.

Covalent and noncovalent strategies for targeting Lys102 in HIV-1 reverse transcriptase.

Reverse transcriptase (RT) is one of three key proteins responsible for the replication cycle of HIV-1 in the host. Several classes of inhibitors have been developed to target the enzyme, with non-nucleoside reverse transcriptase inhibitors forming first-line treatment. Previously, covalent RT inhibitors have been identified and found to bind irreversibly to commonly mutated residues such as Y181C. In this work we aim to circumvent the issue of NNRTI resistance through targeting K102, which has not yet been identified to confer drug resistance. As reported here, 34 compounds were synthesized and characterized biochemically and structurally with wild-type (WT) HIV-1 RT. Two of these inhibitors demonstrate covalent inhibition as evidenced by protein crystallography, enzyme kinetics, mass spectrometry, and antiviral potency in HIV-1 infected human T-cell assays.

Human immunodeficiency virus 1 5'-leader mutations in plasma viruses before and after the development of reverse transcriptase inhibitor-resistance mutations.

Human immunodeficiency virus 1 (HIV-1) reverse transcriptase (RT) initiation depends on interaction between viral 5'-leader RNA, RT and host tRNA3Lys. Therefore, we sought to identify co-evolutionary changes between the 5'-leader and RT in viruses developing RT-inhibitor resistance mutations. We sequenced 5'-leader positions 37-356 of paired plasma virus samples from 29 individuals developing the nucleoside RT inhibitor (NRTI)-resistance mutation M184V, 19 developing a non-nucleoside RT inhibitor (NNRTI)-resistance mutation and 32 untreated controls. 5'-Leader variants were defined as positions where ≥20 % of next-generation sequencing (NGS) reads differed from the HXB2 sequence. Emergent mutations were defined as nucleotides undergoing a ≥4-fold change in proportion between baseline and follow-up. Mixtures were defined as positions containing ≥2 nucleotides each present in ≥20 % of NGS reads. Among 80 baseline sequences, 87 positions (27.2 %) contained a variant; 52 contained a mixture. Position 201 was the only position more likely to develop a mutation in the M184V (9/29 vs 0/32; P=0.0006) or NNRTI-resistance (4/19 vs 0/32; P=0.02; Fisher's exact test) groups than the control group. Mixtures at positions 200 and 201 occurred in 45.0 and 28.8 %, respectively, of baseline samples. Because of the high proportion of mixtures at these positions, we analysed 5'-leader mixture frequencies in two additional datasets: five publications reporting 294 dideoxyterminator clonal GenBank sequences from 42 individuals and six National Center for Biotechnology Information (NCBI) BioProjects reporting NGS datasets from 295 individuals. These analyses demonstrated position 200 and 201 mixtures at proportions similar to those in our samples and at frequencies several times higher than at all other 5'-leader positions. Although we did not convincingly document co-evolutionary changes between RT and 5'-leader sequences, we identified a novel phenomenon, wherein positions 200 and 201 immediately downstream of the HIV-1 primer binding site exhibited an extraordinarily high likelihood of containing a nucleotide mixture. Possible explanations for the high mixture rates are that these positions are particularly error-prone or provide a viral fitness advantage.